Global nuclear structure effects of tensor interaction 
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A direct fit of the isoscalar spin-orbit (SO) and both isoscalar and isovector tensor coupUng 
constants to the /5/2 — /7/2 SO splittings in *'^Ca, ^®Ni, and **Ca nuclei requires a drastic reduction 
of the isoscalar SO strength and strong attractive tensor couphng constants. The aim of this work 
is to address further consequences of these strong attractive tensor and weak SO fields on binding 
energies, nuclear deformability, and high-spin states. In particular, we show that contribution to 
the nuclear binding energy due to the tensor field shows generic magic structure with the tensorial 
magic numbers at N{Z)=14:, 32, 56, or 90 corresponding to the maximum spin-asymmetries in ldc,/2, 
1/7/2 ©2p3/2, I59/2 ©2(^5/2 and l/iii/2 ©2/7/2 single-particle configurations and that these numbers 
are smeared out by pairing correlations and deformation effects. We also examine the consequences 
of strong attractive tensor fields and weak SO interaction on nuclear stability at the drip lines, in 
particular close to the tensorial doubly magic nuclei and discuss the possibility of an entirely new 
tensor-force driven deformation effect. 

PACS numbers: 21.10.Hw, 21.10.Pc, 21.60.Cs, 21.60.Jz 
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I. INTRODUCTION 

The primary goal of energy density functional (EDF) 
methods is to describe ground-state energies of fermion 
systems, i.e., in nuclear physics applications - masses of 
nuclei. The existence of an universal functional describ- 
ing exactly masses of odd, odd-odd and even-even nu- 
clei is warranted by the Hohenberg-Kohn [l[ and Kohn- 
Sham (HKS) theorems. The HKS theorems, however, 
provide no universal rules underlying the construction of 
such functional. The complexity of the functional and 
our lack of knowledge with respect to the in-medium 
strong interaction that governs the structure of finite nu- 
clei make the situation even more difficult. It does not 
permit the determination for any ab initio constraints of 
the nuclear EDF, except for dilute neutron systems [1]. 
It forces us to use effective functionals with coupling con- 
stants fitted directly to the data. Hence, a proper selec- 
tion of empirical data to be used in the process of con- 
straining parameters of the functional becomes, irrespec- 
tive of the form of functional, the key issue for the overall 
good performance of the nuclear density functional the- 
ory (DPT) [4]. 

The typical strategy used to construct the nuclear EDF 
is to start with either the finite-range Gogny [5j or zero- 
range Skyrme [6j effective interactions and construct non- 
local or local functional, respectively, by averaging the 
interaction within the Hartree-Pock (HP) method. The 
datasets used to adjust the free parameters of the the- 
ory are dominated by bulk nuclear matter data and by 
nuclear binding energies of selected doubly magic nuclei 
with much lesser attention paid to the single-particle en- 
ergies (SPE). 

The major reason is the effective mass scaling of the 
single-particle (SP) level density, g, in the vicinity of 
the Fermi energy, ep- In homogeneous nuclear mat- 
ter the SP level density scales according to the simple 



rule: giep) ~* -^gisp)- In finite nuclei the situation is 
slightly more intricate mostly due to r-dependence of the 
TO* (r) . Several authors 0, H, Q analyzed the SP level 
density scaling arguing that the physical density of SP 
levels around the Fermi energy can be restituted only af- 
ter the inclusion of particle- vibration coupling, i.e., by 
going beyond mean- field (MP). This viewpoint is diffi- 
cult to reconcile with the effective EDF theories. In- 
deed, these theories should warrant a proper value of 
the effective mass through the fit to empirical data and 
readjust other coupling constants to this value of to*. 
This should lead to fairly to* independent predictions, 
provided that the (spherical) SPE are calculated from 
the differences between the binding energies in even-even 
doubly-magic cores and the lowest SP states in odd- 
A single-particle/hole neighbors. Within the EDF ap- 
proach, the mean-field or, more precisely, Kohn-Sham SP 
energies computed in e-e doubly magic core serve only as 
auxiliary quantities. 

A new strategy for fitting the spin-orbit (SO) and ten- 
sor parts of the nuclear EDF was recently suggested 
by our group 10]. Instead of performing large-scale 
fits to binding energies we proposed a simple and in- 
tuitive three-step procedure that can be used to fit the 
isoscalar strength of the SO interaction as well as the 
isoscalar and isovector strengths of the tensor interac- 
tion. The entire idea is based on the observation that the 
f-j/2 ~ fb/2 SO splittings in spin-saturated isoscalar '*°Ca, 
spin-unsaturated isoscalar ^^Ni, and spin-unsaturated 
isovector *^Ca form distinct pattern that can neither be 
understood nor reproduced based solely on the conven- 
tional SO interaction. Following the general philosophy 
of the nuclear DPT we compute the f-j/2 ~ h/2 SO split- 
tings from the differences between the binding energies 
of these doubly-magic cores and their odd-^ neighbors. 
However, we use the same functional to calculate both 
the ground state energies, what is well justified, as well 
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as the low-lying SP excitations in the odd-^ neighbors. 
How reasonable is the latter assumption it remains to be 
studied. 

The procedure reveals the need for a sizable reduction 
(from ■--^20% up to ~35% depending on the parameteri- 
zation and, in particular, on the value of m*) of the SO 
strength and at the same time, for much stronger tensor 
fields as compared to the commonly used values. The 
new parameterization systematically improves the per- 
formance of the functional with respect to SP properties 
like the SO splittings or the magic gaps but deteriorates 
the binding energies . The aim of the present work is 
to address further consequences of strong attractive ten- 
sor and weak SO fields on binding energies, time-even and 
time-odd polarization effects, and the nuclear deforma- 
bility. The paper is organized as follows. In Sect. ITTl we 
briefiy present the theoretical background of our model. 
In Sect, mil we show that the contribution to the binding 
energy due to the tensor interaction forms a generic pat- 
tern closely resembling that of the shell-correction with 
the tensorial magic numbers shifted up as compared to 
the standard magic numbers toward iV(Z)=14, 32, 56, or 
90. The tensorial magic numbers reflect the maximum 
spin-asymmetry in Idg/a, I/7/2 © 2p3/2, I59/2 ® 2^5/2, 
and I/111/2 ffi 2/7/2 configurations, respectively, in the 
extreme SP scenario at spherical shape. The tensorial 
magic structure is smeared out by configuration mixing 
caused by pairing and deformation effects. In Sect. IIVI 
we demonstrate that one can construct the EDF capa- 
ble to reproduce reasonably well both the SO splittings 
as well as binding energies of the doubly magic spherical 
nuclei. In Sect. |V] we discuss the infiuence of strong ten- 
sor fields on time-even and time-odd polarization effects 
in the /7/2 — /5/2 SO splittings. In Sect.|Vl]we analyze 
deformation properties of the new functionals and dis- 
cuss a possible novel mechanisms related to the onset of 
nuclear deformation in the presence of strong attractive 
tensor fields. The paper is concluded in Sect. IVIII This 
analysis complements our preliminary results communi- 
cated in two earlier conference publications [ll|, [l^ . 



II. THEORY: FROM TWO-BODY SPIN-ORBIT 

AND TENSOR INTERACTIONS TO ENERGY 
DENSITY FUNCTIONALS AND MEAN-FIELDS 

In our study we explore the local energy density func- 
tional 7i(r). It is the sum of the kinetic energy and the 
isoscalar (t = 0) and isovector {t = 1) potential energy 
terms: 



densities, currents and their derivatives: 



Zm ^ — ' 

t=o,i ^ 



(1) 



which are conventionally decomposed into parts build of 
bilinear forms of either only time-even or only time-odd 



CJptTt + C/J? + C^'ptV ■ 3 



(2) 



(3) 



For the time-even, pt, Tt, and Jt, and time-odd, S(, Tt, 
and jt, local densities we follow the convention introduced 
in Ref. [l^ (see also Refs. [3, [IE] and references cited 
therein) . 

In the present study we focus on the spin-orbit and 
tensor parts of the EDF: 

= C^]5l + C^f,, (4) 
T^^^o = Co^-^poV- Jo + C7VV- Ji. (5) 

These two parts of the EDF are strongly tied together 
through their mutual and unique contributions to the 
one-body spin-orbit potential. This relation can be best 
visualized by decomposing the spin-current tensor den- 
sity J^i/ into scalar J^^^ , vector and symmetric-tensor 
(2) 

densities J^iJ : 



-J^°'>6 



ie 7 



7(2) 



(6) 



and going to the spherical-symmetry (SS) limit where 

the scalar and the symmetric-tensor densities j'i^J 
vanish identically. In this limit the spin-current tensor 
density: 

„ SS limit 1 J tt2 SS limit i -,2 /on 

^ > ^ef_,^riJri and r > (8) 



reduces, therefore, to the spin-orbit vector density having 
a single radial component, Jt = ^Jt(r). The variation of 
the tensor and SO parts of the EDF over the radial SO 
densities J(r) gives the spherical isoscalar {t = 0) and 
isovector (t = 1) SO MFs, 



Wr - ^(cfMr)-C7''^)L-S 



2r 



dr 



(9) 



which can be easily translated into the neutron {q n) 
and proton {q = p) SO MFs, 



W, 



so 



4r 



(10) 



dr 



dr 



Below we perform calculations without assuming neither 
spherical nor time-reversal symmetries. General expres- 
sions for the SO mean-fields can be found in numerous 
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references, see e.g. Refs. [T1,[T1|, and will not be repeated 
here. We provide here spherical formulas (O and pU]) 
because they are crucial for understanding of our fitting 
strategy of the tensor and spin-orbit coupling constants 
which relies on the /7/2 — /5/2 SO splittings in spherical 
doubly magic ^"Ca, ^^Ni and "SCa nuclei, see Ref. [l^ 
and discussion below. 

The functional of the form ([l])-® can be obtained by 
averaging the conventional Skyrme effective interaction, 
vsk{f), within the Skyrme-Hartree-Fock (SHF) approxi- 
mation [T3| . In such an approach twenty EDF coupling 
constants Ct are uniquely expressed by means of ten aux- 
iliary Skyrme- force parameters, see Hence, the SHF 
approximation superimposes relatively strong limitations 
to the EDF. It serves, though, as a reasonable starting 
point for further studies. The approach to the EDF start- 
ing from the effective force indicates also possible ways 
of generalization of the nuclear EDF. In particular, in- 
vestigation of the tensor component requires the use of 
a generalized Skyrme force augmented by a local tensor 
interaction vxir): 

v{r) = vsk{r) + VT{r) (11) 

where 

VT{r) = ^ie| [3(ai • k')(a2 • k') - (ai ■ a2)k'^]S{r) 

+,5(r)[3(^i.k)(<T2-k)-(ai-a2)k2]| 
+io[3(ai • k')5(r)(a2 • k) - (ai • <T2)k' • S{r)k\ ,(12) 

and, conventionally, r = ri — r2 and k — — ^(Vi — V2) 
are relative coordinate and momentum while k' is com- 
plex conjugation of k acting on the left-hand side. By 
averaging vxi^) within the HE approach, one obtains 
the following contribution to the time-even part of the 
EDF 

snr-" = fsf ( Jf V - |sf J.' + hBl E(-^*%)' (13) 

flu 

where 

B^^-lite + 3to) Bl = \{t,-to). (14) 

Note, that there are two independent contributions to the 
tensor part of the EDF. The Skyrme force contributes, 
through the exchange term, to the tensor part of the 
EDF in an uniform manner, i.e., it depends on a unique 
coupling constant (|4]). In contrast, the tensor force con- 
tributes to the EDF in a non-uniform way. Hence, the 
tensor force generates a clear theoretical need to gener- 
alize the EDF (HI) by using three independent coupling 
constants multiplying each of the three terms appearing 
in Eq. Q: 

Cr ( + Ci3,- + Cr E(41)'- (15) 



The effects of such an extension, where the new coupling 
constants need to be adjusted can only be probed in de- 
formed nuclei. Since conventional effective interactions 
and functionals are rather successful in describing nu- 
clear deformation there is no first-hand motivation for 
such a generalization. Hence, in the present study, we 
do not implement this possible extension of the EDF, 
and we use the unique tensor coupling constants C/, as 
defined in Eq. 

The contribution to the time-odd part of the EDF com- 
ing from the tensor interaction is: 

Snf'' = BUsfTt-SsfFt) 

+ Bf^(sfAst + 3(V-St)2), (16) 

where 

Bt^M^e~Q and Sf^ = -^(3ie + to). (17) 

New terms which appear in the time-odd part of the 
EDF, namely the ^ St ■ Ft and ^ (V ■ St)^ will not be 
considered here mostly due to lack of clear experimental 
indicators allowing to fit their strength. Extensive discus- 
sion linking the Skyrme forces to the tensor component 
in the EDF including, in particular, the definition of the 
density Ft can be found in Ref. [lg| , see also Ref. [31 ■ 

The starting point of our consideration is always the 
conventional Skyrme-force-inspired functional with cou- 
pling constants fixed at the values characteristic for either 
SkP [13, SLy4 [13, or SkO [19^ Skyrme parametriza- 
tion. The variants of the EDF with tensor and spin-orbit 
strengths modified along the prescription of Refs. [l^, US] 
will be marked by an additional subscript T: SkPy, 
SLy4T, and SkOr. In the time-odd sector we test two 
variants of the functional with the coupling constants fit- 
ted to the empirical values of the s— wave Landau param- 
eters dll, m, go=OA, (7q=1.2 and to the Gogny-force 
values of the p— wave Landau parameters (2^ . l23j | gi=- 
0.19, gl^O.62: 

go = No{2Q + 2C^PpI'^), g^ = -2NoC^ (3pl'\ (18) 

5(5 = iVo(2Cr + 2CJPpI'^), g[ = -2NoCj Ppl'\ (19) 

where (3 = {'iir'^ /2f/'^ , and N^'^ = -k'^K^ l2m*kF is an 
effective-mass-dependent normalization factor. In these 
variants of the EDF we additionally assume density in- 
dependence of Cf coupling constant, set the spin-surface 
term C^^^^O to zero, and assume gauge-invariant rela- 
tions C{ — —CI and C^^ — C^^ ■ Concerning the time- 
odd tensor coupling constants, , the following two pos- 
sibilities will be tested: (i) broken gauge-symmetry sce- 
nario with Cj fitted to the Landau parameters and (n) 
gauge-invariant scenario with — —C( determined us- 
ing the time-even coupling constants C/. The variants 
of the EDF with spin fields defined using the Landau pa- 
rameters will be labeled either by a subscript £5 or by 
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a subscript Lb for the gauge-invariant and the gauge- 
symmetry-violating functionals, respectively. 



III. TOPOLOGY OF TENSOR CONTRIBUTION 
TO THE NUCLEAR BINDING ENERGY 

Recent revival of interest in the tensor interaction was 
triggered by empirical discoveries of strong and system- 
atic changes in the shell structure of neutron-rich oxy- 
gen [2I] , neon [25I , sodium [H, [23] , magnesium [111 , or 
titanium pol . l3Cll| nuclei including new shell gap open- 
ing at 7V=32. These empirical discoveries were success- 
fully interpreted within the nuclear shell-model after in- 
troducing the so-called monopole shifts. To account for 
the data the monopole shifts are: (i) attractive between 
j| and orbitals and {ii) repulsive between j| and 
orbitals where — I ± 1/2. The physical origin of 
these monopole shifts was attributed to the tensor inter- 
action [31I, m, [11] . Soon after successful shell- model cal- 
culations, the mechanism was confirmed to work within 
self-consistent mean-field calculations using finite range 
Gogny force augmented by the finite-range tensor inter- 
action, see [3J]. It was shown in Ref. [3^ that, apart from 
explaining the shell-structure evolution in light exotic nu- 
clei, the tensor interaction was also capable to account 
for empirical trends in the relative positions of the I/5/2 
and lp3 /2 levels in copper [H, [IB] isotopes or for the evo- 
lution of the single-particle I/111/2 — I37/2 level splittings 
versus N in antimony ^] isotopes. 

The local tensor interaction within the SHF approxi- 
mation was studied first in Ref. [4l|. Based on the SPE 
analysis the effective functional coupling constants C/ 
were evaluated to lie within the triangle Jknown as the 
Brink-Stancu-Flocard (BSF) triangle] marked schemat- 
ically in Fig. [1] For strictly pragmatic reasons, like 
technical complexity and lack of firm experimental in- 
dicators constraining further the BSF estimate, in many 
Skyrme parameterizations (including SIII 43] , SLy4 [l^ , 
SkM* [ii], SkO di]), the tensor terms are simply dis- 
regarded by setting C/ = 0. Moreover, the C/ cou- 
pling constants established through fits to, predomi- 
nantly, bulk nuclear data seem to contradict the BSF esti- 
mates. Indeed, the isoscalar tensor coup ling constants of 
such popular forces hke SLy5 [11], SkP [l3] or Skxc 
are relatively weak while their isovector coupling con- 
stants are positive. These coupling constants lie outside 
the BSF triangle as shown in Fig. [TJ 

On the other hand, direct fits to the SP level split- 
tings [H, [H, [4^] clearly reveal that drastic changes in the 
commonly accepted tensor coupling constants are needed 
in order to accommodate for the SP data. This is vi- 
sualized in Fig. [T] where the new parameterizations are 
marked by black triangles and black diamonds, respec- 
tively. The fact that new values of C/ are still slightly 
scattered is a consequence of different fitting strategies 
as well as different starting point parameterizations used 
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FIG. 1: The isovector, Cf , versus the isoscalar, Cq , tensor 
coupling constant. Open circles mark values representative 
for several popular parameterizations fitted predominantly to 
the binding energies of spherical nuclei. Black diamonds rep- 
resent coupling constants deduced recently from direct fits 
to the SPE and the SP splittings [lol ]. Black triangles rep- 
resent fits of Refs. "3?, "Sy, ^401 . The shaded area shows the 
so-called BSF triangle refiecting the range of the tensorial pa- 
rameters deduced in a pioneering paper by Brink, Stancu and 
Flocard [HI. 



by different groups. 

Our strategy of fitting the coupling constants of the 
nuclear EDF, see Ref. [10], differs from the strategies ap- 
plied by other groups. Unlike the other groups, we fit si- 
multaneously the isoscalar spin-orbit, Cg^"', as well as the 
isoscalar, Cg , and the isovector, C/, tensor couphng con- 
stants using a simple three-step method. The entire idea 
of our procedure is based on the observation that the em- 
pirical 1/7/2-1/5/2 SO splittings in ^^Ca, '^^Ni, and ^®Ca 
form a distinct pattern, which cannot be reproduced by 
using solely the conventional SO interaction. The read- 
justment of the coupling constants goes as follows: (i) 
in the first step Cj^'' is established in the isoscalar spin- 
saturated nucleus *°Ca; [ii) next, the Cq coupling con- 
stant is readjusted in spin-unsaturated isoscalar nucleus 
^^Ni; [iii) finally, the C/ coupling constant is readjusted 
to the spin-unsaturated isovector nucleus '^^Ca. 

Our results, see Refs. [13, [HI [Hj [13] i show that dras- 
tic changes in the isoscalar SO strength and the tensor 
coupling constants are required as compared to the com- 
monly accepted values. In turn, one obtains systematic 
improvements for such single-particle properties like SO 
splittings and magic-gap energies. It is also interesting 
to note that the isoscalar SO and the isoscalar tensor 
coupling constants resulting from such a fit are, to large 
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FIG. 2: (Color online) The isovector (top), the isoscalar (middle) and the total (bottom) tensor contribution to the nuclear 
binding energy obtained from spherical (left) and deformed (right) HFB calculations. Both sets of the calculations were done 
using the SLy4T interaction in the particle-hole channel and the volume-5 interaction in the particle-particle channel. Vertical 
and horizontal lines indicate the single-particle tensorial magic numbers at spherical shape. See text for further details. 



extent, independent on the parameterization and equal to 
C^^ K, -60 ± lOMeVfm^ and « -40 ± lOMeVfm^, 
respectively. The uncertainties are rough estimates re- 
flecting the sensitivity of the method. The isovector ten- 
sor couphng constant, C/, is less certain. It depends on 
the actual ratio Cq'-' /CY'^ of the SO couphng constants 
which, in the adjustment process, was kept fixed to its 
Skyrme force value. This is due to lack of empirical data 
in **Ni which does not allow for firm independent read- 
justment of the fourth coupling constant C^^ ■ 

The infiuence of the tensor interaction on nuclear 
SPE and SP level splittings has been analyzed by many 
authors. It was shown that the tensor interaction 
leaves unique and robust fingerprints when the SPE 
and SP level splittings are studied along isotopic or 
isotonic chains. It appears also that the contribution 
to the binding energy coming from the tensor interac- 



tion, 5Bt {N, Z) , shows several highly interesting and ro- 
bust topological features. In particular, the contribution 
6Bt{N, Z) shows a generic pattern closely resembling 
that of a shell-correction. The tensorial magic numbers 
at iV(Z)=14, 32, 56, or 90 correspond to the maximum 
spin-asymmetries in 1^5/2, l/7/2®2p3/2, lff9/2ffi2d5/2 and 
l/iii/2©2/7/2 single-particle configurations, respectively, 
in the extreme SP scenario at spherical shape. The ro- 
bustness, i.e., model independence of the tensorial magic 
pattern, results from rather unambiguously established 
order of single-particle levels which is relatively well re- 
produced by state-of-the-art nuclear MF models, in par- 
ticular, in light and medium-mass nuclei. Note, that the 
tensorial magic numbers are only slightly shifted as com- 
pared to classic magic numbers at N{Z)=8, 20, 28, 50, 
and 82. 

The tensorial magic pattern is clearly visible in Fig. [21 
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FIG. 3: (Color online) Tensor contribution to the total binding energy obtained from spherical HFB calculations. The map is 
overlaid with the map of known, according to Ref. 4^], nuclei. The figure illustrates these regions where strong tensor effects 
may be expected in neutron- and proton-rich nuclei. 



Panels a)-c) show the contribution to the binding energy 
coming from: a) the isovector part of the tensor term, 
SEf = Ci J d^rJl{r), b) the isoscalar part of the ten- 
sor term, 6Ej — Cq J d'^rjg(r), and c) the total tensor 
contributions to the EDF. These calculations were per- 
formed using spherical Hartree-Fock-Bogolyubov (HFB) 
code HFBRAD |46|] with the SLy4T functional of Ref. [Tg 
in the particle-hole channel and the volume (5-interaction 
in the pairing channel. This part of Fig. [5] shows several 
interesting features including: 

• Additional smearing of the SP tensorial magic 
structure due to configuration mixing caused by 
nuclear pairing. In light and medium-mass nuclei 
a substantial tensor contributions located in rela- 
tively broad regions centered around the SP ten- 
sorial magic numbers. In heavier nuclei, where 
pairing effects are relatively stronger due to larger 
density of SP levels, the erosion of the SP tenso- 
rial magic structure is stronger. The maximum of 
the tensor contribution is shifted away from the SP 
magic numbers. The details, however, are strongly 
model-dependent mostly due to large uncertainties 
in the positions of the SP levels in heavier nuclei. 
Indeed, Skyrme models have persisting problems 
to reproduce absolute positions of the experimen- 
tal SPE as shown recently in Ref. [47]. 

• Contribution due to the isovector part of the tensor 
interaction is much weaker than the isoscalar con- 
tribution. This conclusion depends on the tensorial 



coupling constants. 

• The isoscalar tensor interaction creates oscillatory 
effects in nuclear masses that depend on the de- 
gree of spin-unsaturation in a given nucleus. This 
additional non-uniform N and Z dependence may, 
in particular, obscure conclusions deduced from 
widely used binding-energy indicators technique. 

The second major source of configuration mixing is 
due to the spontaneous breaking of spherical symmetry 
inherent to the MF method. The influence of nuclear 
deformation on the topology of the tensor contributions 
to the binding energy is illustrated in Fig. [5Ji,e,f. The 
calculations presented in the figure were performed for 
e-e nuclei with N > Z ranging from 6 < Z < 64 us- 
ing the HFODD code The same SLy4T interaction 
of Ref. [lOj was used in the p-h channel and volume-(5 
interaction in the p-p channel. It is clearly visible that 
the effect of deformation does not change the topology 
of the tensor energy contribution but strongly reduces 
its magnitude. One should stress though that quantita- 
tive estimate of the deformation effect is uncertain. The 
magnitude of the deformation is extremely sensitive to 
the balance between SO and tensor strengths. This ef- 
fect will be discussed in detail in Sect. fVl 

Fig. [3] shows again the total contribution to the bind- 
ing energy calculated using the spherical HFB model. 
The map is overlaid with the map of known, according 
to Ref. [i^, nuclei. The aim of the figure is to illus- 



7 



^ 3 
> 

W 

u 



1 I 1 I 1 

SkpT SLy4 

p.... 


1 1 1 1 


: oM — 1 

E>0 * 

1 I 1 1 1 


L \ ^ 



"°«cr''»N?'-..z7^".4n'''' 
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binding energies in spfierical doubly magic nuclei and in ^''Zr. 
The calculations have been done using the SkOy' (white dots) 
and SkOy'/ (gray dots) functionals. The conventional SLy4 
Skyrme-force result (black dots) is also shown for the sake of 
comparison. 



trate mass-regions where enhanced tensor effects, and 
in turn perhaps new physics, may be expected on the 
neutron-rich and proton-rich side. On the neutron-rich 
side the regions of interest, i.e., those which are or can be 
accessible experimentally in the nearest future, include: 
Z « 14 and TV « 32, Z « 32 and « 56, and Z w 56, 
N « 90. In particular, recent measurements of exotic 
f2Mg28 and f§Al29 by Bauman et al. [i^ (see also the 
discussion in Ref. [50j) approach closely the first of the 
above mentioned mass regions. However, whether or not 
extra binding due to the strong attractive tensor interac- 
tion gives rise to stabilization of these nuclei and nuclei 
around them remains to be studied. Mean-field calcu- 
lations using conventional Skyrme forces predict these 
nuclei to be bound i51j. 



IV. ENERGY DENSITY FUNCTIONAL FITTED 
TO THE SINGLE-PARTICLE SPIN-ORBIT 
SPLITTINGS AND TO THE TOTAL BINDING 
ENERGIES OF SPHERICAL NUCLEI 

The topology of the tensor contribution to the total 
binding energy is, as discussed above, a generic feature re- 
lated to shell-structure and the degree of spin-saturation. 
The quantitative features including the total magnitude 
and the isovector to isoscalar ratio of the tensor contri- 
butions depend, however, upon the actual values of the 
tensorial coupling constants. The two strategies of fitting 
effective forces, namely the conventional one based on the 
large-scale fit to the binding energies and the one based 
on the fit of C/ and Cj'' coupling constants directly 
to the SO splittings seem to yield contradicting results. 
This is clearly visible in in Fig. [1] where the ranges of the 
C/ strengths deduced using these two methods simply 



exclude each other. 

A direct fit to the SO splittings leads to drastic changes 
in the isoscalar SO strength and the tensor coupling con- 
stants as compared to the commonly accepted values. In 
particular, the pronounced reduction of the SO strength 
which varies from 20% for large effective mass, m*, forces 
to even 35% for low (m* « 0.7) forces imperils the agree- 
ment with observed binding energies. This effect cannot 
be solely compensated by strong attractive tensor fields. 
Hence, further readjustments of the other coupling con- 
stants of the EDF are necessary to assure good perfor- 
mance on masses. 

In Ref. [l3| we have demonstrated that a considerable 
improvement in the quality of mass fits can be achieved 
be relatively small readjustments of the EDF coupling 
constants. For the purpose of this work we have per- 
formed similar calculations using multi-dimensional min- 
imization technique but starting from the SkPy force. In 
both cases tiny modifications (of the order of a fraction 
of a percent) in the coupling constants clearly improve 
the quality of the mass fit as compared to the SLy4T and 
SkOT forces but still do not provide the quality of the 
original SLy4 and SkO parameterizations. 

Inherent to the multi-dimensional minimization tech- 
nique is a merit function being a subject of minimiza- 
tion. In our calculations the merit function is constructed 
out of relative deviations from measured masses of se- 
lected spherical doubly magic nuclei. It appears that the 
merit function in the multi-dimensional space spanned 
by the EDF coupling constants varying in the minimiza- 
tion process is very steep for some specific directions and 
extremely flat in others. It implies that the entire min- 
imization problem is not well defined. We will visualize 
this by taking as a starting point the SkO parameteriza- 
tion. 

By reducing the SO strength by 15% corresponding 
to Ci « -65.6MeVfm5 and C( w 84.5MeVfm5 and 
taking Ci = -44.1McVfm5 and C( = -91.6MeVfm5 
we create a modified version of the SkOr parameteriza- 
tion of Ref. [l^ . This parameterization, dubbed SkOr' , 
has slightly stronger SO term and more attractive tensor 
fields as compared to the SkOr. These changes aim to 
improve mass performance of the SkOr- Direct calcula- 
tions show, see Fig. m that the SkOr' functional repro- 
duces masses at a similar level of accuracy as one of the 
most popular SLy forces. The calculations illustrated in 
Fig. |4] were performed using the code HFODD with 20 
spherical shells. We were forced to use the HFODD code 
because the spherical HFBRAD code has no two-body 
center-of-mass correction implemented. 

The performance of the SkOr' force can be further 
improved in many different ways. One example, dubbed 
SkOT" , is illustrated in Fig. |31 This force was obtained 
by readjusting isoscalar and isovector central fields in 
the following way: 1.00015C^ and Cf ^ 0.99Cf. 

As a result, standard deviation drops from (T(SkOT') ~ 
1.663 MeV to (T(SkOT") ~ 1.475 MeV where, for compar- 
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ison, (T(SLy4) w 1.879 MeV. Similar improvements can 
be made by, for example, readjusting the density depen- 
dent term. 

A reasonable performance of the SkOy' or SkOy" func- 
tional with respect to the binding energies of spherical 
doubly-magic nuclei is of great interest. It may help to re- 
solve the conflict concerning the tensor and SO coupling 
constants resulting from (local) fits to the SO splittings 
and the SPE [lOjESi preferring strong tensor and weak 
SO terms on one hand and from large-scale fits to the 
binding energies pointing toward weak tensor and 
stronger SO terms on the other hand. Indeed, our present 
result indicates that one should explore in large-scale 
fits functionals having non-standard forms including, in 
particular, functionals having non-conventional isovector 
spin-orbit which characterize the SkO functional. 



V. TIME-EVEN AND TIME-ODD 
POLARIZATION EFFECTS IN THE PRESENCE 
OF STRONG TENSOR FIELDS 

The aim of this section is to analyze the polarization 
phenomena in the presence of strong tensor fields and 
to show that, in spite of the relatively large readjust- 
ments as compared to the SkOr parameterization, the 
SkOr' functional can still reproduce the empirical I/5/2 
and 1/7/2 SO splittings quite well. Fig. [5] illustrates the 
neutron (left hand side) and proton (right hand side) SO 
splittings between the I/5/2 and I/7/2 SO partners in 
^°Ca, ''^Ca, and ^^Ni. The empirical data are marked 
by filled circles. The values shown are average means of 
empirical results taken from Refs. [sl, [s^ (see also Ta- 
ble III in Ref. [13] for compilation of the empirical SO 
splitting data). Error bars represent deviations from the 
mean- values. 

Open circles illustrate the results of our calculations 
using the original SkO parameterization. Open and filled 
diamond represent calculations using SkOr' functional. 
In all variants of the calculations the Cf^" strength was 
set to zero to assure convergence. Different panels rep- 
resent different variants of the calculations concerning 
treatment of the time-odd sector. Panels a) and e) show 
bare, unpolarized SO splittings deduced directly from the 
SP spectra calculated in the doubly magic '^^Ca, ^^Ca, 
and ^^Ni nuclei. Results presented in the other panels 
are calculated from the binding energies in the doubly 
magic cores and their one-particle(hole) odd-j4 neighbors 
following the prescription given in Ref. [lo| . In the cal- 
culations the odd-A binding energies correspond to fully 
aligned {ly) = j states at oblate, for one-particle, and 
prolate, for one-hole, nuclei, respectively. Unlike in our 
previous study the present calculations include po- 
larization effects exerted by the odd particle or hole on 
the even-even core in the presence of strong attractive 
tensor fields. 

In order to visualize the role of the tensor interaction, 
in particular in the time-odd sector, we performed three 



different variants of the calculations. Figures b) and f ) 
include only time-even (mass and deformation) polariza- 
tion effect. These results were obtained by setting all 
time-odd coupling constants to zero. Panels c) and g) 
illustrate the effect of the time-odd fields. These results 
include both the time-even and time-odd polarizations. 
In these calculations we use gauge invariant functional 
with the = — C/ tensor coupling constants fitted to 
the SO splittings. All other coupling constants in this 
run, except C^'^ = 0, are equal to the values given by 
the SkO parameterization. Note that the splittings cal- 
culated in this way match almost perfectly the empirical 
data. Note also, that the time-odd polarization effects 
are indeed large, reaching a few hundred keV. 

Figures d) and h) show two variants of the calcula- 
tions with spin fields readjusted to match empirical Lan- 
dau parameters, namely the SkOr'Ls and the SkOT'L^ 
(see Sect. ini). The SkOr'Lg variant is labeled by open 
diamonds. It corresponds to fully gauge invariant func- 
tional. In this variant we use the = — C/ tensor cou- 
pling constants fitted to the SO splittings and the spin- 
field coupling constants readjusted to the s-wave Landau 
parameters according to Eqs. p^ - P^ . 

In the SkOT'Ls variant, which is labeled by filled di- 
amonds, we readjust first the Cj^ coupling constants 
to the Gogny values of the p-wave Landau parameters 
gi ~ —0.19 and g'l = 1.2. This leads to gauge-symmetry 
violating functional with the tensorial coupling constant 
C5' « -44.1 MeV fm^ and « -91.6 MeV fm^ in the 
time-even channel and Cj ~ 9.2 MeV fm^ and « 
— 29.9MeVfm^ in the time-odd channel which are used 
subsequently to calculate spin-fields coupling constants. 
Note, that the SO splittings are quite sensitive to the way 
the functional is set up in the time-odd sector. It means 
that the entire concept of fitting the time-odd coupling 
constants to the Landau parameters in the presence of 
the strong tensor terms must be reconsidered. In particu- 
lar, the p-wave parameters deduced from the Gogny force 
gi — —0.19 and g'l — 1.2 lead to the coupling con- 
stants which are completely inconsistent with the time- 
even values C/ deduced from the SO splittings. In turn, 
the SkOT'Ls and SkOT'L^ have entirely different spin 
fields with coupling constants equal Cg w 426.4 MeV fm^ 
and Cf « 48.6 MeV fm^ and w 18.0 MeV fm^ and 
Cf w 155.9 MeV fm^, respectively. 



VI. EFFECT OF TENSOR FIELD ON 
NUCLEAR DEFORMATION 

In the previous section we have shown that one can 
construct the EDF capable to reproduce reasonably well 
binding energies of the spherical doubly magic nuclei 
and, simultaneously, account for the SO splittings of the 
1/7/2 — 1/5/2 SP levels. By inspecting closer the results 
presented in Fig. 2 one observes very strong time-even 
(mass and deformation) polarization effect on the calcu- 
lated SO splitting in ^^Ni - compare Fig. [5^ and Fig. [2h 
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FIG. 5: SO splittings between the neutron and vlfy/i (left part) and the proton and iTlf-j/2 (right part) orbitals 

in *''Ca, **Ca, and ^''Ni nuclei. Filled (open) circles mark empirical and theoretical splittings calculated using the conventional 
SkO force, respectively. Open and filled diamonds indicate calculations performed using the SkOy/ functional. Panels a) and 
e) show bare, unpolarized SO splittings deduced directly from the SP spectra calculated in the doubly magic ''"Ca, **Ca, and 
^^Ni nuclei. Parts b) and f) include only time-even (mass and deformation) polarization effect. Panels c) and g) include both 
the time-even and time-odd polarization effects. Finally, figures d) and h) show two variants of the calculations with spin fields 
readjusted to match empirical Landau parameters, namely the SkOy'L^ (open diamonds) and the SkOj-'i^ (filled diamonds) 
functionals. 



or Fig. [2^ and Fig.[2t. The effect is definitely stronger for 
the SkOr' than for the SkO, indicating that the reduction 
of the spin-orbit combined with strong attractive tensor 
fields can affect deformation properties which are rather 
well captured by conventional Skyrme forces. Hence, nu- 
clear deformability in the presence of strong attractive 
tensor fields and reduced SO potential poses very strin- 
gent test for this new class of functionals. 

The aim of this section is to show that deformation 
properties indeed depend sensitively on the balance be- 
tween the SO and tensor fields. We will discuss two con- 
trasting cases. First, we will consider an example of the 
yrast superdeformed (SD) band in ^^Ni [sj]. The band is 
formed by promoting two protons and two neutrons from 
the 1/7/2 to the I/5/2 orbital or, in terms of more appro- 
priate from mean-field point of view asymptotic Nilsson 
model quantum numbers, from [303]7/2 to [321]3/2 Nils- 



son orbital as illustrated schematically in Fig. [6^. 

In the spherical ground state of ^^Ni the entire I/7/2 
is fully occupied while the I/5/2 is empty. This creates 
large spin-asymmetry and, in turn, large contribution to 
the ground state due to the tensor field. By promot- 
ing four-particles from the I/7/2 to the I/5/2 orbit one 
creates the SD state which has reduced spin-asymmetry 
as compared to the ground state. The reduced tensor 
field shifts the SD state up in energy with respect to the 
ground state. The subsequent reduction of the SO po- 
tential shifts the I/5/2 orbit and, in turn, the entire SD 
band down in energy. This compensating mechanism is 
illustrated in Fig. [5}D,c,d. The figures show the results 
of the self-consistent quadrupole-constrained HF calcu- 
lations for the ground state (GS) and the SD configura- 
tions. Since it is impossible to go diabatically through the 
GS-SD configuration crossing region the self-consistent 
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FIG. 6: Part a) shows schematically the mechanism underlying the formation of the SD band in ^''Ni. It is associated with 
4p-4h isoscalar excitation from the I/7/2 to the I/5/2 orbital. Figures b)-d) show potential energy curves for the GS and the 
SD bands calculated using the SkO, SkOTX, and SkOy' functionals, respectively. Parts b) and c) demonstrate the effect of 
strong attractive isoscalar tensor field and part d) the effect of reduced spin-orbit field on the excitation energy of the SD band. 
Part e) shows the change in the tensor energy associated with the 4p-4h excitation leading to the SD band. 



results for the GS and SD configurations were, for the 
sake of simplicity and clarity, extrapolated diabatically 
through this region. This does not affect the physics dis- 
cussed below. 

Fig.O) shows the calculations performed using the con- 
ventional SkO force. These calculations predict the 0+ 
SD state to be excited by ~4 MeV with respect to the GS 
what agrees quite well with the empirical estimate 
Readjustment of the tensor coupling constant in the SkO 
to the value characteristic for the SkOy/ functional (this 
functional is called SkOrx ) shifts the position of the 0+ 
SD state by 4 MeV up in energy. This intermediate step 
is illustrated in Fig.[Sl;. The change in the tensor field on 
the passage from the GS to the SD minimum is shown 
in Fig. [6^. Finally, Fig. [6ji shows the result obtained 
using full SkOr/ functional. The effects of the reduced 
SO and strong attractive tensor fields almost cancel each 
other restoring the position of the 0+ SD state close to 
its empirical (and close to the SkO) value. 

The second example is shown in Fig. [71 The figure il- 
lustrates potential energy curves in *°Zr calculated using 
the HFB model with volume-delta interaction, VS{r), of 
the strength V"=-190MeV in the particle-particle chan- 
nel. The three curves represent the SkO, SkOyx, and 
SkOT' functionals in the particle-hole channel. All curves 



are normalized to the spherical minimum in order to fa- 
cilitate further discussion. 

Unlike the ^^Ni, the ^°Zr is spin-saturated in the spher- 
ical minimum. Building up deformation is associated in 
this nucleus with partial occupation of the 139/2 sub- 
shell. It leads to increasing spin asymmetry and, in turn, 
to extra attraction due to the tensor terms. This effect is 
clearly visible for the SkOrx functional, see Fig. [T] The 
mechanism is strong enough to create a well deformed 
minimum. In this case, reduction of the SO strength 
shifts the 1^9/2 sub-shell up with respect to the negative 
parity fp levels. In turn, the well deformed minimum 
is also lifted up in energy ending up slightly higher as 
compared to the SkO prediction. 

These two examples show that the SkOy' functional 
have deformation properties quite similar to the conven- 
tional Skyrme functionals at least in isoscalar N Ki Z nu- 
clei. The situation is slightly more intricate in the isovec- 
tor channel due to the uncertainties of the C^"' strength 
and, in turn, in the C( coupling constant. Nevertheless, 
the two cases analyzed above clearly show that: 

• Nuclear deformation properties strongly depend on 
the balance between tensor and spin-orbit terms 

• Detailed and systematic studies of nuclear defor- 
mation in the presence of strong tensor fields open 
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FIG. 7: Potential energy curves calculated using quadrupole- 
constrained Skyrme HFB model. Solid curve represents the 
SkO calculations. Dotted and dashed curves illustrate the 
SkOTX and SkOy' results, respectively. Note, that tensor 
field tends to favor strongly elongated shapes leading to well 
deformed absolute minimum. The subsequent reduction of 
the SO strength, however, shifts the entire potential energy 
curve up in energy close to its original SkO position. 

up new, promising venue which may help to tune 
the efltective tensor coupling constants C/. 

• There is a large, and so far completely unexplored, 
potential to modify deformation properties by ex- 
tending the tensor term from the uniform form used 
here to the non-uniform form (jlSp. 

• Possible modifications of collective rotational mo- 
tion, which is inherently related to the spontaneous 
breaking of spherical symmetry within the EDF for- 
malism, opens up yet another almost completely 
unexplored area in studying tensor fields. 

VII. SUMMARY AND CONCLUSIONS 

The direct fit of the isoscalar spin-orbit and both 
isoscalar and isovector tensor coupling constants to the 
/5/2 - /7/2 SO splittings in ^OQa, ^^Ni, and ^^Ca requires 

(i) a drastic reduction of the isoscalar SO strength and 

(ii) strong attractive tensor coupling constants [1^. In 
this work we address the global nuclear structure conse- 
quences of this novel fitting strategy of the nuclear EDF. 
Among others, we show that contribution to the nu- 
clear binding energy due to the tensor field shows generic 
magic structure with the tensorial magic numbers at 
N{Z)— 14:, 32, 56, or 90 corresponding to the maximum 
spin-asymmetries in 1^5/2, I/7/2 ® 2^3/2, I59/2 © 2^/2 
and I/111/2 ©2/7/2 single-particle configurations and that 
these numbers are smeared out by pairing correlations 
and deformation effects. 



We explicitly construct the functional, dubbed SkOr', 
which is able to reproduce simultaneously the /2 ~ /? /2 
SO splittings in ^°Ca, ^^Ni, and ''^Ca nuclei and the bind- 
ing energies of spherical nuclei. In fact, one can con- 
struct many parameterizations reproducing these data 
in a more or less equivalent manner. This is due to 
the fact that multi-dimensional merit-functions which are 
minimized in the fitting process are flat in certain direc- 
tions and the entire minimization procedure is not well- 
defined. Reasonable performance on the binding ener- 
gies of the SkOr' functional which is characterized by its 
non-conventional isovector SO term is very interesting 
by itself. Indeed, this result may open up a possibility 
to bridge the and C^'^ coupling constants resulting 
from (local) fits to the SPE with the values resulting from 
global large-scale fits to the binding energies by exploring 
non-standard local functionals. 

Using the SkOr' functional we analyze polarization ef- 
fects exerted by one particle and one hole on the values 
of the /5/2 — /7/2 SO splittings. We show that the po- 
larization effects are large and very sensitive to the way 
the functional is set up in the time-odd channel. Fits 
to the Landau parameters are uncertain due to rather 
poorly known p-wave Landau parameters. In particular, 
the use of the p-wave Landau parameters deduced from 
the Gogny force, which is advocated in Ref. [2^, leads 
to strong gauge-symmetry violating effects and, in turn, 
large differences in the coupling constants between 
the gauge-invariant SkOr' Ls functional and the gauge- 
symmetry violating functional SkOT'L^. 

We also demonstrate that deformation properties in 
atomic nuclei can be easily and strongly modified in the 
presence of strong tensor fields and that these proper- 
ties are extremely sensitive to the balance between the 
tensor and the SO coupling constants. We show that, 
in the particular case of the SkO^' functional, the ten- 
sor effects are almost perfectly compensated, at least in 
the isoscalar channel, by the reduced SO potential. We 
suggest that the role of the tensor interaction, in par- 
ticular in the time-odd channel, can be studied through 
dynamical effects induced by fast nuclear rotation. For 
this purpose one needs to select and use nuclear states 
representing, as close as possible, an unperturbed single- 
particle motion in order to suppress other effects or cor- 
relations which may obscure the analysis. The examples 
include superdeformed bands which are known to be very 
well described u sing simple one-body cranking approxi- 
mation, see Ref. [55j and references therein, or terminat- 
ing states [20j . 
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